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ABSTRACT 

A  multimedia  model  is  developed  and  applied  to  selected  organic  chemicals 
in  evaluative  and  real  environments.   The  model  employs  the  fugacity 
concept  and  treats  four  bulk  phases:  air,  water,  soil,  and  bottom  sediment, 
which  consist  of  eubcompartments  of  varying  proportions  of  air,  water,  and 
nineral  and  organic  matter.    Chemical  equilibrium  is  assumed  to  apply 
within  (bat  not  between)  each  bulk  phase.   Expressions  are  included  for 
emissions,  advective  flows,  degrading  reactions,  and  inter-phase  transport 
by  diffusive  and  non-diffusive  processes.    The  solution  of  the  steady- 
state,  constant  emission  condition  is  a  relatively  simple  algebraic 
expression.    Input  to  the  model  consists  of  a  description  of  the 
environment,  the  physical-chemical  and  reaction  properties  of  the  chemical, 
end  the  emission  rate.  The  model  is  applied  to  six  chemicals  in  the  region 
of  Southern  Ontario  and  the  calculated  concentrations  compared  with 
observations.  The  results  suggest  that  the  model  may  be  used  to  predict 
the  likely  environmental  fate  of  chemicals  and  estimate  concentrations  for 
defined  locations  and  at  defined  times  within  an  order  of  magnitude. 


Sommaire 

On  elabore  un  modele  portant  sur  plusieurs  phenomenes  pour 
l'appliquer  a  certaines  substances  organiques  en  milieux 
experimental  et  reel.   Ce  modele  fait  appel  au  concept  de  la 
fugacite  et  porte  sur  quatre  phases  globales  de  la  matiere  : 
1'air,  l'eau,  le  sol  et  les  sediments.   Chacune  de  ces  phases 
est  composee  de  sous-categories  selon  sa  teneur  en  air,  en  eau, 
et  en  substances  minerales  et  en  matieres  organiques.   On 
presume  1' existence  d'un  equilibre  chimique  a  1'interieur  de 
chacune  des  phases  globales,  mais  non  entre  celles-ci.   Le 
modele  comprend  des  formules  pour  les  emanations,  l'advection, 
les  effets  de  degradation  et  le  transfert  de  masse  par 
diffusion  et  par  d'autres  procedes.   La  solution  pour  les 
emanations  en  regime  constant  est  presentee  sous  forme 
d' expression  algebrique  relativement  simple.   Le  modele  fait 
1' analyse  des  donnees  sur  1 'environnement,  sur  les  proprietes 
physiques  et  chimiques  de  la  substance  et  sur  le  debit 
d' emission.   On  applique  le  modele  a  six  composes  chimiques  que 
l'on  retrouve  dans  le  sud  de  1' Ontario  pour  ensuite  comparer 
les  resultats  obtenus  aux  concentrations  observees.   Les 
resultats  semblent  indiquer  que  le  modele  permet  de  prevoir  le 
sort  probable  des  substances  chimiques  dans  1' environnement  et 
d'estimer,  selon  un  ordre  de  grandeur,  leur  concentration  a  des 
endroits  et  a  des  moments  precis. 


INTRODUCTION 

Environmental  models  of  varying  degrees  of  sophistication  have  been 
developed  to  describe  or  predict  the  multimedia  fate  and  behavior  of 
organic  contaminants  (1.2.8.4).  These  models  maybe  steady- or  unsteady- 
state  In  nature  and  employ  the  chemical's  physical  chemical  properties, 
reactivity,  transport  characteristics,  and  extent  of  release  to  the 
environment  to  produce  a  comprehensive  picture  of  the  chemical's  pathways 
between,  and  ultimate  concentrations  in,  various  media,  such  as  soil,  air, 
water,  and  sediment.   This  modelling  capability  is  a  necessary  adjunct  to 
toxicity  testing  in  assessing  risk  from  chemical  exposure   and  for 
determining  priorities  for  research  and  regulatory  action.  It  can  be 
applied  to  new  and  existing  chonlcals. 

Of  concern  are  chemicals  that  are,  or  will  be.  in  frequent  commercial  use 
and  therefore  can  be  expected  to  be  discharged  constantly  to  the 
environment  over  a  period  of  time.  Examples  are  pesticides,  PCB's.  wood 
preservatives,  and  by-products  of  incineration  and  dry-cleaning  operations. 

We  develop  here  a  novel,  simplified,  fugacity-based  model  to  predict  this 
steady-state  behaviour.  The  model  treats  four  bulk  phases:  air,  water, 
soil  and  sediment,  which  consist  of  subcompartments  of  fluid  (air  and/or 
water),  particulate  (aerosols,  suspended  particles,  and  organic  matter)  and 
other  phases  (fish).  Equilibrium  partitioning  of  chemicals  is  assumed  to 
apply  within  each  of  the  four  bulk 


Expressions  for  environmental  emissions  and  advection  .  into,  reaction  and 
advection  from,  and  diffusive  and  non-diffusive  transfer  between  these  bulk 
phases  are  incorporated  in  a  set  of  mass  balance  equations  as  previously 
described  (2).  Reducing  the  number  of  key  environmental  media  to  four 
permits  a  single  algebraic  solution.  Unlike  previous  models,  non-diffusive 
transfer  processes,  such  as  wet  and  dry  atmospheric  particulate  deposition, 
rain  washout,  leaching  to  groundwater  from  soil,  and  sediment  resuspension, 
are  included. 

Pugacity 

In  previous  models,  partitioning  of  mass  between  phases  has  been  described 
by  the  equilibrium  criterion  of  fugacity  as  an  alternative  to  chemical 
potential.  Fugacity  can  be  considered  to  be  the  partial  pressure  of  a 
chemical  in  a  phase  and  is  related  to  its  concentration  in  that  phase  by 
the  expression 

C  =  fZ  ID 

where  C  is  concentration  (mol/aPj,  f  is  fugacity  (Pa)  and  Z  is  the  fugacity 
capacity  (mol/m3.Pa).  Wien  equilibrium  exists  between  two  phases,  the 
chemical  fugacities  are  equal  and  partitioning  can  be  described  by  the 
ratio  of  their  Z  values,  that  is 


K12  =  C!/C2  -  fZi/fZ2  -  Zi/Z2 


(2) 


where  K12  Is  the  dimensianless  equilibrium  pertiticn  coefficient.  The 
calculation  of  Z  values  for  environmental  phases  has  been  reviewed  by 
Mackay  and  Paterson  (1) . 

Briefly,  Z  for  air  is  calculated  using  the  gas  law,  and  Z  values  for  other 
phases  are  calculated  from  it  using  estimated  partition  coefficients  such 
as  sorption  partition  coefficients,  Henry's  Law  constants,  and 
bioconcentration  factors.  Partitioning  of  a  chemical  between  the  organic 
phases  of  soil,  sediment,  or  fish  and  water  can  be  described  in  terms  of 
its  octanol/water  partition  coefficient  (K^)  and  the  fraction  organic 
content  in  that  phase.  For  soil  or  sediment,  the  partition  coefficient  is 
0.41  Kc^x  L/kg  where  x  is  the  organic  content  of  the  phase  (1).  Fish  to 
water  partitioning  is  described  by  the  relationship  0.048Kow  (8)« 

Chemicals  move  between  phases  by  diffusive  and  non-diffusive  processes. 
Diffusive  transfer  rates  between  two  phases,  1  and  2,  can  be  described, 
using  the  fugacity  equilibrium  concept  as 

N  -  D12fi  -  I>2if2  <3> 

where  N  is  the  flux  between  the  two  phases  (mol/h) ,  D12  and  D21  are 
transfer  coefficients  with  units  of  mol/h.Pa,  and  fj  and  f2  are  the  phase 
fugacities.  The  difference  between  fj  and  f2  determines  the  direction  of 
diffusive  flux  which  takes  place  from  high  to  low  fugacity.  The 
coefficients  D12  and  D2i  are  functions  of  Z  values,  interfacial  areas  and 


diffusion  properties  In  adjacent  phases  (2),    In  the  case  of  simple, 
reversible  diffusion,  Dj2  and  D21  are  equal  and 

N  -  D12(fi  -  f2)  (4) 

It  has  been  shorn  (2)  that  D  may  be  expressed  as  the  product  of  kAZ  where  k 
Is  a  mass  transfer  coefficient  and  A  Is  the  interfacial  area.  In  some 
cases,  there  may  be  two  diffusive  resistances  in  series,  each  characterized 
by  an  individual  mass  transfer  coefficient.  D  may  also  be  expressed  as 
BAZ/Y,  where  B  Is  a  diffusivity  and  Y  Is  the  effective  diffusion  path 
length,  B/Y  being  equivalent  to  k  (3) . 

Jfon-diffusive  or  one-way  transfer  processes  between  phases,  for  example, 
wet  or  dry  particle  deposition  from  the  atmosphere  to  soil  or  water,  or 
suspended  sediment  deposition  or  resuspension,  can  also  be  described  by  a 
transport  paramenter  as 

N  ■  GC  »  GZf  -  Df 

where  G  is  volumetric  flow  rate  (ro3/h) ,  and  D  again  has  units  of  mol/Fa.h. 

Similarly,  advective  out-flows  in  air  and  water  can  be  described  by  a  flew 
rate  (G  n^/h) ,  concentration  product. 

rate  of  advection  ■  GC  =  GZf  =  DAf 
or  rate  of  advection  *  k^VZf  ■  D^f 


ttoere  1^  is  a  first  order  rate  constant  equal  to  G/V.  Advectlve  Inflow  can 
be  Included  as  0%  or  DAfB,  in  *#iich  subscript  B  refers  to  the  background 
inflow  concentration. 

First  order  reaction  processes  in  a  phase  can  also  be  described  by  a  D 
value  as 

rate  of  reaction  «*  kRVC  «  kgVZf  =  Dpf 

where  1^  Is  a  first  order  rate  constant  (fcT*l,  V  is  the  phase  volume  (m3). 
and  Dr  again  has  units  of  nol/Pa.h. 


MODEL  DESCRIPTION 

1.  Environmental  Media 

Four  major  bulk  compartments  are  defined  each  of  total  volune,  VA,  -   air, 

water,  soil,  and  sediment  -  subscripted  1-4.   These  phases  consist  of  a 

combination  of  subcoopartments  of  pure  and  particle  phases  of  volune,  Vjj, 

and  are  defined  as  a  volune  fraction,  $|j  or  Vij/Vj,  of  the  bulk  phase. 

For  example,  bulk  water  consists  of  "pure"  water,  suspended  solids,  and 

fish;  soil  consists  of  mineral  matter,  organic  matter,  water,  and  air;  air 

consists  of  pure  air  and  aerosols;  sediment  is  mineral  matter,  organic 

matter,  and  water. 

It  is  assumed  that  equilibrium  exists  within  each  bulk  phases,  ie.  a  common 
fugacity  applies.  This  Implies  that  the  times  required  to  reach 
equilibrium  within  a  bulk  phase  are  short  compared  to  times  required  to 
reach  equilibrium  between  bulk  phases. 

This  environment  can  be  evaluative  in  nature  as  described  in  previous 
papers  (3,6)  or  can  represent  an  actual  region.  We  apply  the  model  here  to 
i)  an  environment  scaled  to  represent  Southern  Ontario  and  11)  an 
evaluative  environment  of  area  1  sq  km  with  similar  phase  proportions.  The 
air  height  has  been  reduced  to  2000  m  from  6000  m  to  reflect  conditions  of 
atmospheric  accessibility  In  Southern  Ontario.  The  compartment  volumes  for 
the  two  scenarios  are  given  in  Table  1. 


Table  1 

3  3  2 

Properties  of  the  Phases:  Volume  Fractions,  Organic  Carbon  (OC)  Fractions,  Densities  (kg/m  ),  Volumes,  (m  ),  Areas  (m  ), 

and  Depths  (m) 


Bulk  Phase 
(J) 


Air 
(1) 


Subcompartments  (j) 
Volume  Fractions 
Water    Solids    Biota 
(2)       (3)      (4) 


Fraction   Density 
OC  in 
solids 


Volumes  Areas  Depths 

S.  Ontario   Evaluative  S.-. Ontario  Evaluative 


\lr    (1) 

1.0 

0 

5x10" 

•12 

0 

— 

1.19 

Cater    (2) 

0 

1.0 

5x10" 

6 

lxlO-6 

0.2 

1000 

Soil    (?) 

0.2 

0.3 

0.5 

0 

'0.02 

1500 

iediment    (4) 

0 

0.7 

0.3 

0 

0.04 

1420 

)ensity 

1.19 

1000 

2400 

1000 

4x10 


4x10 


1.2x10 


14 


12 


10 


8x10 


8 


2x10' 

2xl07 

4 
6x10 

4000 


2x10 


8x10 


11 


10 


12x10 


8x10 


10 


10 


1x10 


4x10" 


6x10" 


4xl0_ 


2xl0J 
50 
0.1 
0.01 


00 


For  each  subcompartaent,  a  Z  value  or  fugaclty  capacity  Is  calculated. 
These  values  are  then  combined  with  the  volume  fractions  to  give  an  overall 
Z  value  for  the  bulk  phase  as  illustrated  in  Table  2,  ie.  ZA  is  the  sun  of 
the  ^ijZij  terns. 

I)  Bulk  Air 

Vapor  and  particulate  natter  (aerosols)  are  assumed  to  be  the  principal 
components  of  the  bulk  air  phase.  From  the  gas  law,  the  fugaclty  capacity 
for  vapour  or  pure  air,  Zu,  is  calculated  as  1/RT  where  R  is  the  gas 
constant  (8.314  Pa.mS/mol.K)  and  T  is  the  absolute  temperature.  Particles 
in  air  occupy  a  volume  fraction.  fa  of  2X10"11,  equivalent  to  about  50 
ug/m3.  Partitioning  between  particulate  and  gaseous  air  phases  is 
inversely  proportional  to  the  vapor  pressure  (7)  with  the  particulate-air 
partition  coefficient  expi-essed  as 

K13  *  6x106/Pls  -  Z13/Zn 

where  pSL  is  the  sub-cooled  liquid  vapor  pressure  (Pa) .    Therefore  the 
fugaclty  capacity  for  particulates  in  air  is  given  by 

z13  »  exioe.Zn/pSL 

and  the  overall  Z  value  is  calculated  as  shown  in  Table  2. 


»- 


11)  Balk  Hater 

Hater  Is  considered  to  consist  of  a  pure  water  phase  containing  suspended 
sediment  particles  and  fish.  The  fugacity  capacity  for  pure  water  Z22  is 
defined  as  the  reciprocal  of  the  Henry's  Law  constant  or  1/H  (1).  The  Z 
values  for  the  particle  (Z23)  and  biota  or  fish  (Z24)  phases  are  calculated 
»  PKij./>ij-Z22  (or  K^Z^)  where  PKi;J  is  the  particle  (fish) /water 
partition  coefficient  (L/kg).   Kjj   is   the   dimensionless  partition 
coefficient,  and  f> ^    is  density  of  the  particulate  or  fish  phases  (kg/L) 
(1).  In  the  calculation  of  PK23  a  value  of  0.2  is  assumed  for  mass 
fraction  organic  content  (x1;j).    The  fish/water  partition  coefficient, 
PK24,  is  0.048  Kow-   These  Z  values  are  combined  with  their  appropriate 
volume  fractions,  ^23  and  ^24  of  5x10-6  (equivalent  to  about  12  ng/L)  and 
10-6,  respectively,  to  give  an  overall  fugacity  capacity  as  showi  in  Table 
2.   it  should  be  noted  that  the  concentration  of  "fish"  is  unreal  is  tlcally 
high  and  represents  the  inclusion  of  other  biota,  especially  micro- 
organisms. 

lil)  Bulk  Soil 

Following  Jury  et  al  (8) ,  soil  is  assumed  to  contain  a  soil  particle  phase 
(volume  fraction  /33  of  0.5)  combined  with  air  (^31  of  0.2)  and  water  (jfe2 
of  0.3) .  The  fugacity  capacity  for  the  solid  phase,  Z33,  is  calculated  as 
described  earlier  from  the  product  of  the  pure  soil/water  partition 
coefficient  and  the  Z  value  for  water  as  K32.Z22.  The  fraction  organic 
content  in  the  solid  phase  x33  is  0.02,  and  the  density  ^33  is  2.4  kg/L. 
The  Z  values  of  soil,  air,  and  water  are  combined,  as  shorn  in  Table  2, 
with  their  volume  fractions  to  give  an  overall  fugacity  capacity. 


(I 


Table  2 

Definition  of  Z  values 


Subcompartments  {j) 

Air  Hi) 
Water  ( ±2  J 


1/RT         R 
1/H  or  Cs/Ps  H 

cs 

PS 


Soil,  Sediment  Solids,   XijKoc  ij/H   Xj_j 
&  Water  Particles  U3) 

Koc 


Reference 


Aerosols  {13) 


Bulk  Phases  (i) 

Air  (1} 
Water  (2) 
Soil  (3} 
Sediment  (4) 


pa 


=  8.314 (Pa  m3/rool  K)  T=absolute  temp  K 

m  Henry's  Law  constant  (Pa  m3/mol) 

=  aqueous  solubility  (mol/m3) 

=  vapor  pressure  (Pa) 

=  fraction  organic  carbon 

■  organic  carbon  partition  coefficient 
=  0.41  Kqw 

=  density  (kg/L) 


6x106/PlSRT 


z:  =  Zu  +  &3Z13 
Z2  =  Z22  *  023z23  +  024z24 
z3  ■  ^31Z11  +  #32z22  +  ^33z33 
Z4  =  ^42"22  +  #43z43 


(1) 
(1) 


(5) 


(7) 


Biota  (fish)    (i4)  0 . 048  oi4Kbw/H        Kqh  =  octanol  water  partition  coefficient  (6) 


lv)  Bulk  Sediment 

The  bulk  phase  is  assaoed  to  consist  of  two  subccnpartnenta,  pure  sediment 
arri  sediment  pore  water  with  volume  fractions  of  ^43  and  ^42,  with  values 
of  0.3  and  0.7,  respectively.  As  in  the  case  of  soil,  the  2  value  for  pure 
sediment  (Z43)  is  calculated  as  K42.Z22.  where  X43  -  °'°*  and^  ■  2-* 
kg/L,  and  the  combined  bulk  Z  value  is  deduced  as  in  Table  2. 

Transfer  Parameters 

Assumed  transfer  parameters  necessary  for  calculation  of  intermedia  fluxes 
are  given  in  Table  3,  areas  being  given  in  Table  1.  These  parameters  are 
converted  into  D  values  as  shown  in  Table  4.  A  summary  of  their 
derivations  follows. 

i)  Air-Water 

Chemicals  move  between  air  and  water  by  diffusive  and  non-diffusive 
processes.  Diffusive  processes  of  volatilization  and  absorption  are 
described  by  the  parameter  Dy  where 

l/Dy  *  1/Dva  +  1/DyW  "  l/t^VA^H )  +  1/(^^12222) 

where  kyA  and  kw  «**  air-  and  water-side  KTC  (mass  transfer  coefficents) 
(m/h) ,  Dva  and  Dvw  tbe  transfer  parameters  (mol/h.Pa) ,  and  A12  is  the 
interfacial  area  (m*).  The  reciprocals  of  these  D  terms  can  be  considered 
resistances  which  are  in  series  and  combine  to  give  an  overall  air-water 
transfer  resistance  1/Dy.  The  same  mass  transfer  coefficients  are  assumed 
to  apply  to  all  chemicals  with  values  of  3  m/h  and  0.03  m/h  for  kVA  and 
kw,  respectively.    These  values  are  slightly  reduced  from  those  used 


Table  3 

Assumed  Transport  Parameters 


n 


Parameter 


Symbol 


Value 


Air  side  MTC  over  water  kVA 

Water  side  MTC  kVW 

Transfer  rate  to  higher  altitude  Us 

Rain  rate  (m3rain/m2area.h)  Uq 

Scavenging  ratio  Q 

Dry  deposition  velocity  Up 

Air  side  WTC  over  soil  ^5A 

Diffusion  path  length  in  soil  Y3 

Molecular  dif fusivity  in  air  B^a 

Molecular  diffusivity  in  water  %i 

Water  runoff  rate  from  soil  Uww 

Solids  runoff  rate  from  soil  Ug^ 

Water  side  MTC  over  sediment  kyw 

Diffusion  path  length  in  sediment  Y4 

Sediment  deposition  rate  UdX 

Sediment  resuspension  rate  UrX 

Sediment  burial  rate  Ugx 

Leaching  rate  from  soil  ^L 

to  ground  water 


1  m/h 

.01  m/h 

.01  m/h  (90  m/y) 

9.7X10"5  m/h  (.85m/y) 

200,000 

10.8  m/h  (.003  m/s) 

1  m/h 

.05  m 

.04  m2/h 

4.0xl0-6  m2/h 

*3.fcxl0-5  m/h  ( ."54m/Y) 

2.3X10-8  m/h  (0.0002  m/y) 

.01  m/h 

.005  m 

4.5xl0"8  m3/m2.h(.00O4  m/y) 

l.lxlO-8  m3/m2.h( .0001  m/y) 

3.4xl0-8  m3/m2.h( .0003  m/y) 

3.9xlO-5  m3/m2.h(.34  m/y) 


t  *1 


Table  4 

Calculation  of  D  Parameters 
Coapartfflent      Process 


Individual  D 


Total  D 


alr(l)  -  water(2) 


diffusion  Dv    -  l/(l/kVAAl2Zll  +  l/kwAl2Z22> 

rain  Dqh  ■  Z22  •  *k& 

wet  deposition  Dew  ■  °DH   •  Q  •  ^13  ■  z13 

dry  deposition  Dpw  -  Gpw  •  ^13  •  z13 


D12<vHiqw+DDw«>Pw 

D2ia€>v 


air(l)  -  soll(3)        diffusion 


Dg    =  l/(l/kSAA13Zn  + 

Y3/(Ai3(BA3Zn  +BM3Z22)) 
rain  Dqs  =  GqS  *  Z22 

wet  deposition  Dqs  "  ^DS  •  Q  •  F1&13 
dry  deposition  Dps  =  Gps  •  013   ■  z13 


D13-DS+DQS+ODS+DPS 
D31<>S 


soil{3)  -  water(2) 


soil  runoff        Dsw  =  Gjc  •  z33 
water  runoff      D^  ■  Gw  .  Z22 


D32"DgWfDWW 
D23-0 


sediment  (4)-water(2)dif  fusion  Dy    =  l/U/kjmta^  ♦  Y4/W4A24Z22) 

deposition         B&  -  <%«  •  ^3  D24'«Y+DDX 

resuspension      Drx  =  Grx  .  Z43  D42"DY+°RX 


reaction  either  bulk  phase      Drj  =  fcfeiViZi 

or  pure  phase      Drj-j-  kRij   .  VijZAj 


kRiViZi 
£(kR13VijZij] 


advection 


bulk  phase  Dp^  =  GjZi 


GiZi 


previously  to  account  for  a  larger  average  chemical  molecular  weight  and  to 
correspond  more  closely  to  values  reported  for  lakes. 


The  non-diffusive  processes  of  chemical  transfer  by  rain  dissolution  and 
wet  and  dry  atmosphere  deposition  to  water  are  described  by  the  appropriate 
Qd2i  products  as  D(j,,  DDw.  and  Dpw,  respectively.  The  method  of 
calculating  of  these  D  values  is  given  in  detail  by  Mackay  et  al  (3). 
Briefly,  rain  which  falls  at  a  rate  Uq  of  9.7xl(T6  (mSrain/n^area.h)  or  m/h 
or  0.85  m/year  on  water  contains  dissolved  chemical  in  equilibrium  with  the 
air  and  with  a  concentration  of  Cq  mol/ta3.  The  flux  from  air  to  water  in 
the  dissolved  phase  Kgw  is  expressed  as 

Xq»  =  UQfc12CQ  -  G^.Cq 
"  Dgw-fl  (mol/h) 

where  Gqw  is  the  volumetric  flow  rate  of  rain  to  the  water  surface  (m3/h) . 
E^j  is  thus  Gqw.Z22,  since  Cq  is  Z22*l- 

For  wet  particulate  deposition,  each  rain  drop  is  considered  to  scavenge 
all  particles  from  an  atmospheric  volume  Q  (200,000)  times  larger  than 
itself.  Since  the  particles  occupy  a  volume  fraction  $^3.  the  volumetric 
flux  on  water  of  chemical  in  association  with  wet  particles,  %*  becomes 

Hpw  »  Uq.Q.^i3.A12.Cp 
■  %W-Cp 
c  °DW-fl 


t* 


where  Dr**  is  G^.Z^,  %jw  is  the  wet  deposition  rate  of  particles  on  the 
water  {m3/h} ,  and  CP  is  the  chemical  concentration  on  the  particles. 

For  dry  deposition  of  particles  on  water  the  flux  is 

Npw  =  Up.jta.A^.Cp 

-  Gpw-Ziafi 

Where  Gpw  (m3/h)  is  the  product  of  Up,  the  dry  particle  deposition  velocity 
(assumed  to  be  10.8  m/h  or  0.003  m/s},  the  particle  volume  fraction,  and 
the  water  area. 

The  transfer  parameters  necessary  to  calculate  the  net  flux  between  bulk 
air  and  water  become: 

air-water  D12  =  Dy  +  Dqw  +  °DW  +  DPW 

water-air  D2i  =  Dy 

with  total  flux  N12  =  E>i2fI  -  D2if2 

ii)  Air-Soil 

The  bulk  soil  is  viewed  as  a  matrix  10  cm  deep  of  air,  water,  organic  and 

mineral  matter.  Diffusion  from  soils  can  take  place  in  the  air  or  water 

phases  or  in  parallel  in  both  and  is  usually  described  by  Fick's  first  law 


:>B^'A/  Ail  A. 
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where  B  is  an  effective  dlffusivity  In  air  (Bai)  or  water  (B^i)  (m2/h),AC 
is  the  concentration  difference  (mol/m3)  and  AY  is  the  diffusion  path 
length  (m)  to  which  the  difference  applies.  B  is  calculated  as  in  the 
model  developed  by  Jury  et  al.  (8)  using  the  Miliington  and  Quirk 
expressions  and  the  appropriate  molecular  dlffusivity  and  phase  volume 
fraction.  For  soil  air,  the  effective  diffusivity  is 

where  B^  is  air  molecular  diffusivity,  /31  is  volume  fraction  of  soil  air 
and  (1-^33)  is  the  porosity  of  the  soil. 

Similarly,  for  diffusion  in  soil  water  (**iich  is  loportant  for  soluble,  low 
volatility  chemicals),  the  effective  dlffusivity  is 

The  D  values  for  the  parallel  diffusive  flows  in  soil  air  and  water  are 

DAD  *  BA3.Ai3.Zn/Y3   air 
Efcn  =  Bw3.Ai3.Z22/Y3   water 

BA3/Y3  a"3  Bvo/Ys  can  be  considered  equivalent  to  mass  transfer 
coefficients  in  soil  air  and  water,  Y3  being  the  soil  diffusion  depth  of  5 
cm  or  half  the  total  depth. 


m  the  Interests  of  simplicity,  cannon  molecular  diffusivities  acre  used  for 
all  chemicals  of  0.04  in  air  and  ixlO-6  (mVh)  in  water. 

A  mass  transfer  coefficient  kg^  (of  1  m/h)  combined  with  the  air  fugacity 
capacity  and  area  characterises  the  air-side  diffusion  parameter  and 

Recalling  that  the  resistance  l/Dg^  is  in  series  Kith  the  two  parallel 
resistances,  l/D^n  and  l/Dfc  total  D  for  soil-air  diffusion  becomes 

D&  =  1/{1/DSA  ♦  i/tDAD  +  °M))) 

Non-diffusive  transfer  of  chemical  by  atmospheric  deposition  is  included  to 
give  total  D  values  of 

air-soil  D13«Ds  +  I>DS  +  I)PS  +  DQS 
soil-air  D31  ■  D5 

where  D^,  Dps,  and  Dqg  are  parameters  for  wet  and  dry  deposition  and  rain 
wash  out  fran  air  to  soil  and  are  calculated  as  for  the  air-water  exchange 
but  adjusted  to  reflect  the  soil-air  interfacial  area. 


[*l, 


iv)  Water-Sediment 

Vfeter-sediment  diffusion  is  calculated  similarly  to  soil-air  diffusion. 

The  effective  diffusivity  in  sediment  pore  water  is  given  by 

Bw4  -  Bw  .  ^423*33/C1  "  ^43)2 
=  B^  .  /421'33 

and  diffusion  in  pore  water  is  represented  by 

DXW  =  BW4  •  A24  ■  Z22/Y4 

where  Y4  the  effective  depth  is  5  mm. 

The  total  diffusive  transport  coefficient,  Dy,  between  sediment  and  water 
is 

r^  =  i/(i/pxw  +  i/Dyw) 

where  diffusion  in  the  water  boundary  layer  is  described  by  the  parameter 

Dyw  =  KyW.A24-222 
and  %g  is  the  water-side  mass  transfer  coefficient  with  a  value  of  .0.1 


The  one  way  processes  of  transfer  by  suspended  sediment  deposition  and 
resuspenslon  are  described  by  "velocities"  Unx  {4.5x10-3)  and  Ifcx  (l.lxKT 
8)  a3  solids/m2  area.h  or  n/h  which  are  combined  with  sediment  area  to  give 
flow  rates  Gnx  and  Oft*.  «*  respective  transfer  parameters  become 

Djk  *  (%k.Z23 
Djqj  m  QRX-343 

Inclusion  of  these  parameters  gives  overall  D  values  for  sediment-water 
transfer 

water-sediment   D24  »  Dy  +  Dnx 
sediment-water   D42  *  Dy  +  Drx 

The  sediment  depth  has  been  reduced  to  10  mn  which  combines  with  a  burial 
rate  of  .0003  v&/v?  y  of  solids  or  1  ma/year  of  sediment  of  porosity  0.7  to 
give  a  characteristic  burial  time  of  10  years,  which  is  reasonable  for 
freshwater  lakes. 

v)  Soil-Water 

Transfer  of  chemical  from  soil  to  water  is  considered  to  take  place  in  one 
direction  only  and  includes  transfer  by  water  run-off,  Djw.  •"*  associated 
soil  loss  Dgw  calculated  from  appropriate  GZ  values,  ie 

tyw  =  <*l-z22 
DSW  -  GsZ33 


V 


The  parameters  for  total  transfer  become 

soil-water  D32  =  D^w  +  DSW 
water-soil  D23  =  0 

The  water  and  solids  run  off  rates  from  soil  are  expressed  per  unit  area  of 
soil  and  represent  fairly  heavy  run  off  conditions.  The  solids  runoff  rate 
of  .00023  m/y  is  chosen  to  balance  the  sediment  burial  rate  of  .0003  m/y, 
and  the  water  runoff  is  assumed  to  be  40%  of  the  rain  rate  of  0.85  m/y. 
These  values  are  typical  of  precipitation  and  soil  run-off  rates  reported 
by  Environment  Canada  and  the  U.S.  EPA  for  the  Great  lakes  Basin  (9) . 

Reaction  Parameters 

As  stated  previously,  reactive  and  advective  processes  are  characterized  by 

the  parameters  Dr  and  DA  with  units  of  mol/Pa.h. 

The  model  permits  reaction  in  either  the  bulk  or  pure  phase  of  each  media 
but  not  in  both  simultanebusly.   The  reaction  rate  constant  is  combined 

with  the  appropriate  V  and  Z  values.  For  example,  for  reaction  in  bulk  air 


or  for  reaction  in  air  p^-ticles 


**  *<      '        -"-  '    i      ' —  ^  . 
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If  reaction  occurs  In  more  than  one  pure  phase,  the  D  values  are  sunned  to 
give  a  total  for  the  bulk  phases,  that  is 

%Ti  -£»*ijvij2ij 

Advectlan  Parameters 

Mvection  is  into  and  out  of  the  bulk  phases  of  air  and  water.  However,  it 
is  possible  to  incorporate  an  advective  flow  in  particulate  Batter.  Inflow 
rates  and  residence  times  for  air  and  water  are  given  in  Table  5. 

Transfer  to  an  altitude  greater  than  2  km,  leaching  from  soil  and  sediment 
burial  are  included  as  advective  processes.  Transfer  to  higher  altitude 
takes  place  at  approximately  90  m/y  (or  .01  m/h) ,  as  discussed  previously 
(3) .  The  leaching  rate  from  soil,  which  has  been  included  in  an  attempt  to 
detect  chemicals  which  are  potential  groundwater  contaminants,  is  set  at 
40*  of  the  rain  rate  of  0.85  m/y.  Sediment  burial  is  assumed  to  take  place 
at  0.3  mm/year  and  is  equal  to  the  difference  between  suspended  sediment 
deposition  and  resuspension  rates.  The  velocities  are  included  in  Table  3 
and  the  residence  times  in  Table  5. 

We  emphasize  that  these  values  have  been  selected  largely  on  the  basis  of 
being  typical  of  the  Southern  Ontario-Great  Lakes  region.  They  can  easily 
be  modified  for  other  regions. 

The  only  differences  between  the  regional  and  evaluative  models  are  the 
volumes,  areas,  and  advective  flow  rates.  Since  all  other  rates  are 
expressed  on  a  unit  area  or  volume  basis,  they  apply  to  both  models. 


Table  5 

Advection  Rates  end  Residence 


Volume 

Flow  Rate 

Residence 

(m3) 

(n?/h) 

Tine  (h) 

Air 

S.  Ontario 

4X10*4 

3.3x1012 

120 

Evaluative 

2x109 

1.7xl07 

120 

Meter 

S.  Ontario 

4X1012 

3.3x10s 

12000 
(500  days) 

Evaluative 

2x10? 

1700 

12000 
(500  days) 

Transfer 

S.  Ontario 

to 

higher 
altitude 

Evaluative 

Leaching 

S.  Ontario 

from  soil 

Evaluative 

Sediment 

S.  Ontario 

burial 

Evaluative 

2.1x10*     22  years 


lXlO* 

22  years 

4.66X106 

23.3 

2740 

10  years 

0.014 

10  years 

AH 


Emissions 

Emissions  %*  (nol/h)  may  be  Into  any  or  all  of  the  bulk  phases. 

Mass  Balance  Equations 

Mass  balance  equations  are  set  up  for  the  four  bulk  phases  which 
incorporate  emissions,  transfer  between  adjacent  phases,  and  degradatlve 
and  advective  flows.  The  equations  take  the  form: 

%i  +  DAifBi  -£jCfinij  -  fjDji)  +  nAi*i  +  ^Rifi 
or   Nei  +  DAifBi  -  *iljl>ij  +  Dm  +  DRil  +XjtDjifjl  *  ° 

The  solution  is  algebraic  as  shown  in  Figure  I.  Prom  the  bulk  phase 
fugacities,  the  concentrations  and  amounts  for  all  phases  are  calculated  as 
well  as  transfer  and  reaction  rates. 


xt 


Figure  I.  Algebraic  Solution  of  Model  Equations 

f2=X,/Y, 
f,-J,+K,ft 

*4=  J4  +  K4T2 

where    l(  =  D,2  tD,3t  DA!+  DRI+  DST 

|2=  D21  +D24+Da2+Dr2 
l3=D3,  +  D32+DR3+  DL 

J,  =  [NET,  +  D3,J3]/[l|-  D3,K3]       K|  =  D2|/  [l|  -D31K3] 
J3=  NE3/I3  K3=D,3/I3 

J4=  NE4/I4  K4=D24/U 

X,=NET2+  J,D|2+    [J3+JIK3]D32+J4D42 

Y|  =  l2-K|  D|2~  K|K3D32"  K4D42 

NET,  =  NE,+  DA1fBi 

NET2=NE2+  DA2fB2 


ai 


MODEL  OUTPUT 

The  regional  awl  evaluative  model  outputs  Include  the  Input  chemical 
properties,  the  dimensions  and  properties  of  all  bulk  and  subuaiparUnents, 
input  reaction  rate  constants,  input  emission  rates,  and  any  background 
(advective  inflow)  concentrations.  All  fugacitiee  are  provided,  as  are 
concentrations,  in  mol/m3  and  conventional  units.  Transfer  rates  and  D 
values  are  given  for  individual  and  combined  processes.  A  total  mass 
balance  can  be  assembled,  as  shown  f or  KBe  in  Figure  2,  from  **iich  the 
overall  chemical  persistence  is  calculated. 
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%=  percentage  of  total  amount 


Figure  2.   Mass  Balance  for  PCBs  in  the  Great  Lakes 


DISCUSSION 

Models  such  as  this  cannot  be  "validated"  In  the  sane  sense  that  a  simpler 
physical  or  chemical  model  can  be  validated.  Environmental  concentrations 
vary  in  time  and  space,  and  few  data  are  usually  available.   Reaction  rate 
constants  vary  diurnally  and  seasonally.    Emissions  are  rarely  known 
accurately.  There  are  variations  and  uncertainties  in  transport  rate 
parameters,  such  as  deposititon  and  resuspenslon  rates.  Few  environments 
are  at  steady  state  as  described  here.  The  best  that  can  be  hoped  for  is 
that  the  model  corresponds  with  order  of  magnitude  fidelity  to  fragmented 
observations  for  a  variety  of  chemicals  of  quite  different  properties  and 
pathways.   The  chemical-to-chemical  variation  should  be  described  entirely 
by  the  chemicals'  properties  and  emission  rates  and  not  by  adjustable 
parameters  describing  environmental  processes.  But  even  this  modest  model 
capability  Is  invaluable  for  predicting  the  behavior  of  new  chemicals  and 
for  estimating  order  of  magnitude  concentrations. 

As  part  of  the  model  development,  the  model  was  run  for  selected  chemicals 
of  varying  properties,  environmental  concentrations  were  calculated,  and 
compared  with  observations.  A  common  problem  is  estimation  of  emission 
rates.  In  some  cases,  it  is  necessary  *o  back-calculate  what  the  emissions 
may  be  or  must  be  to  create  prevailing  concentrations,  then  check  if  these 
emission  rates  are  in  reasonable  accord  with  rates  of  production  or  use. 

The  chemicals  selected  and  their  properties  are  listed  in  Table  6. 
Emissions,  background  inflow  concentrations .  and  degradation  rate  constants 
are  listed  in  Table  7.  In  the  case  of  benzene,  benzo(a)pyrene  (BAP),  and 
2,3,7,8  TCDD,  production  data  were  used  to  estimate  emissions.  In  the  case 
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Tabis  b.    Ft-yikal  Vesical  PrspsrUes  of  Selected  ChesicsU 
so!,  pt.  %f,  Ec!  y  ioq  Kow 


6-CB  (10) 
tefliBW 

sirex 
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of  6-CB,  hexachlcirobiphenyl  (HCB),  and  alrex,  emissions  were  estimated  on 
the  basis  of  known  uses  and  disposal  methods.  Rate  constants  ware  taken 
from  the  literature  and  judgement  was  used  to  select  a  value  where  a  range 

of  data  existed. 

Table  8  lists  the  calculated  bulk  compartment  concentrations  and  gives 
reported  values  for  comparison.  Table  9  summarises  the  predicted  phase 
concentrations  and  fugacities. 


PCBs 


For  PCBs  (specifically  represented  in  this  case  by  a  hexa-chloro- 
biphenyl),  the  primary  sources  are  water  advection  inflow  (mainly  the 
Niagara  River)  of  1  ng/L  and  air  advection  at  1  ng/B?,  but  with  additional 
local  emissions  into  soil  and  water  due  to  transformer  leakage,  waste 
disposal,  and  seepage  from  landfill  with  a  small  amount  being  released  to 
the  atmosphere  through  combustion. 

The  conditions  shown  are  believed  to  represent  those  prevailing  in  the 
Southern  Ontario  region  In  the  mid-1980s.  Much  higher  concentrations  were 
believed  to  have  been  experienced  in  the  early  1970s.  For  example,  gull 
egg  concentrations  were  then  about  a  factor  of  ten  higher  than  current 
values. 

The  dominant  fate  processes  are  atmospheric  deposition,  volatilization,  and 
sediment  deposition.  Reported  concentrations  of  total  PCBs  are  in  good 
agreement  (ie.,  within  a  factor  of  five)  with  model  estimates,  but  it 
should  be  noted  that  it  is  possible  to  adjust  emission  rates  to  improve  the 
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FCSe  are  a  very  convenient  calibrating  chemical  for  models  because 
they  migrate  widely  throughout  the  environment  and  are  well-studied. 


Production  of  benzene  as  a  chemical  and  fuel  In  Southern  Ontario  Is 
rsported  to  be  approximately  6x10*  kg/year.  It  Is  assumed  that  its  use  as 
an  industrial  chemical  and  fuel  results  in  a  IX  loss  to  the  environment,  of 
which  90*  Is  probably  into  water  due  to  waste  disposal  and  10*  into  air  due 
to  fuel  combustion  and  evaporation.  An  advective  flow  In  air  of  2  ug/m? 


Predicted  values  for  air  and  water  concentrations  are  low  by  a  factor  of 
10.  There  are  no  reported  data  for  soil  and  sediment,  due  probably  to  Its 
low  concentrations  caused  by  its  short  estimated  half-life  of  5  days. 

It  is  probable  that  the  low  concentration  reflects  the  prevalence  of 
measurements  in  urban  areas  where  there  is  a  much  higher  than  average  fuel 
use.  A  good  match  of  model  estimates  with  observations  Is  obtained  if  it 
is  assumed  that  emissions  are  a  factor  of  ten  higher.  Oils  could  represent 
emission  of  say  90*  of  the  benzene  in  10*  of  the  region's  area.  The 
estimated  sediment,  soil,  and  fish  concentrations  are  very  low,  thus  it  is 
unlikely  that  monitoring  these  compartments  will  be  useful. 


Benzo(a)pyrene 

Benzo(a)pyrene  (BAP)  is  mainly  an  atmospheric  pollutant  which  is  deposited 
in  association  with  particulate  matter,  resulting  in  widespread 
contamination  of  soil  and  sediments.  The  emissions  of  65*  into  air  and  35% 
into  soil  assume  additional  direct  soil  contamination  from  aviation  and 
automotive  fuels.  An  air  advective  inflow  concentration  of  1  ng/m3  is 
based  on  data  reported  by  Canviro  (17).  BAP  has  been  the  subject  of 
modelling  by  Ryan  and  Cohen  (18). 

As  with  benzene,  it  is  likely  that  regions  in  which  there  is  a  higher  than 
average  fuel  usage,  especially  diesel  fuel  usage,  experience  concentrations 
which  are  larger  by  a  factor  of  ten  than  those  estimated  by  the  model. 

Hexachlorobenzene 

The  main  route  of  entry  of  hexachlorobenzene  (HCB)  into  the  environment  is 
in  waste  water  streams,  soil  application  of  pesticides,  and  spillage.  It 
is  highly  hydrophobic  and  sorbs  strongly  to  suspended  sediments  with 
subsequent  deposition  to  sediment.  It  is  bioconcentrated  by  fish. 

Emissions  are  assumed  to  be  80%  into  the  water  and  10%  each  into  soil  and 
air.  This  results  in  a  relative  distribution  that  is  in  agreement  with 
reported  concentrations.  The  absolute  value  of  the  emissions  were  chosen 
to  produce  levels  of  the  correct  order  of  magnitude.  It  is  not  known  if 
these  emissions  are  correct. 


Ml  rex 

Mirex  has  been  used  as  a  pesticide  and  a  flame  retardant.  Its  main  routes 
of  entry  to  the  Great  lakes  have  been  via  the  Niagara  and  Oswego  Rivers 
from  production  plants  in  the  U.S.  It  has  been  extensively  monitored  in 
Lake  Ontario.  Sediment  and  suspended  sediment  concentrations  are  well 
quantified,  but  water  concentrations  are  now  generally  below  1  ng/L. 
»teasur«ents  have  been  made  by  Yin  (25)  of  2  pg/L.  iterry  and  Chan  (19) 
estimated  loadings  of  13-20  kg  per  year  to  Lake  Ontario  in  association  with 
suspended  matter.  Good  agreement  was  found  between  predicted  and  reported 
water  concentrations  for  these  emissions.  Sediment  concentrations  are 
underpredicted,  which  is  probably  due  to  previous  contamination. 

Kirex  is  interesting  because  it  represents  a  situation  in  which  there  was  a 
high  emission  rate  in  the  past,  but  present  emissions  are  small.  If  the 
ndoel  is  run  for  two  emission  scenarios  corresponding  to  concentrations  of 
Inflow  water  differing  by  a  factor  of  600,  then  good  agreement  is  obtained 
with  concentrations  experienced  In  the  1970s  and  with  prevailing 
concentrations.  The  steady-state  model  can  thus  be  used  to  estimate 
unsteady  state  behaviour,  provided  that  account  is  taken  of  the  response 
time  of  the  system,  which  in  this  case  is  12.5  years. 


2,3.7,8  Tetrachloro  Dibenzo-p-Dioxir.  (TCDD) 

Polychlorinated  dibenzo-p-dioxins  form  a  goup  of  75  congeners  vtfiich  enter 
the  environment  as  by-products  of  combustion,  waste  disposal,  and  pesticide 
manufacture.  They  are  variable  in  toxicity,  with  2,3,7,8  TCDD  being  the 
most  toxic.  Direct  measurements  of  2,3,7,8  TCDD  are  not  available,  but  it 
is  estimated  that  approximately  40  kg/year  enter  the  Ontario  region  as 
2,3,7,8  TCDD  "toxic  equivalents"  (20).  It  is  assumed  that  80S  of  this 
amount  is  into  the  air  and  10%  each  into  soil  and  water.  It  is  included 
here  primarily  to  illustrate  the  applicability  of  the  model  to  "new" 
chemicals.  Atmospheric  transport  and  deposition  in  association  with 
particulate  matter  are  the  dominant  processes.  The  agreement  between 
estimated  and  the  few  reported  concentrations  is  good. 

Trichloroethylene  (TCE) 

TCE,  a  major  industrial  solvent,  is  introduced  mainly  into  the  air  due  to 
its  high  vapor  pressure.  A  small  percentage  (estimated  here  as  1%)  is 
introduced  to  the  aquatic  environment  through  wastewater  and  disposal. 

The  Ontario  emissions  are  scaled  from  U.S.  production  data  (22).  A 
background  atmospheric  inflow  of  80  ng/m3  }S  included. 

The  resulting  concentrations  are  in  excellent  agreement  with  data  from  La 
Jolla,  California,  and  Liverpool,  England,  as  well  as  the  model  predictions 
of  Cohen  and  Ryan  (22). 


DISCUSSION 

The  model  provides  a  relatively  simple,  rapid  method  of  establishing  the 
relative  importance  of  environmental  fate  processes  for  specific  chemicals 
using  a  minimum  of  input  data.  It  is  believed  to  give  an  adequate 
characterization  of  the  dominant  phases  of  accumulation.  For  example, 
mlrex  is  62*  in  soil  and  35*  in  sediments,  while  benzene  is  46*  in  air  and 
53*  in  water.  These  percentages  are  controlled  by  the  emission  sources  and 
the  chemicals'  properties.  An  overall  persistence  or  residence  time  is 
also  estimated  which  ranges  from  6  days  for  benzene  to  12  years  for  mlrex. 
This  persistence  is  a  function  of  partitioning,  reaction  rates,  and 
advection  rates.  The  relative  importance  of  advective  inflow  and  local 
emissions  can  be  estimated. 

The  evaluative  1  sq  to  version  can  be  used  to  determine  the  dominant 
environmental  partitioning  pathways  and  the  persistence  of  a  chemical, 
using  illustrative  rather  than  real  emissions.   The  computed  absolute 
concentrations  have  no  significance.    Since  many  jurisdictions  have 
interest  in  estimating  the  environmental  fate  of  chemicals  and  the  regions 
corresponding  to  these  jurisdlcations  have  differing  ratios  and  natures  of 
air,  water,  soil,  and  sediments,  it  is  useful  to  have  a  common  evaluative 
test  system  which  can  be  used  by  various  Jurisdictions  to  inter-calibrate 
environmental  assessments.    The  combination  of  evaluative  and  region- 
specific  environments  thus  provides  a  convenient  system  for  international 
assessment  of  chemicals. 


A  useful  feature  of  both  models  is  that  a  sensitivity  analysis  is  possible 
by  which  the  effects  of  varying  properties  of  the  chemicals  and  the 
environment  can  be  explored.  This  establishes  which  properties  should  be 
estimated  most  accurately  and  can  provide  an  estimate  of  likely  error  in 
concentration  as  a  function  of  error  in  property. 

If  emission  estimates  are  available,  the  model  can  be  used  to  estimate 
concentrations  with,  it  is  believed,  order  of  magnitude  (factor  of  10) 
accuracy.   It  is  most  reliable  for  substances,  such  as  PCBs,  which  are 
persistent  and  widely  dispersed  in  the  environment.  The  model  can  be  used 
as  an  interpretive  tool  to  explore  temporal  and  spatial  variations  in 
concentrations.  For  example,  it  appears  that  mirex  has  been  discharged  at 
two  time  periods  at  emission  rates  which  have  differed  by  a  factor  of 
approximately  500.  Since  the  model  is  linear,  concentrations  are  linearly 
related  to  emissions,  thus  there  is  no  need  to  run  the  model  repeatedly  to 
explore  such  changes.  Certain  contaminants,  such  as  benzene  and  BAP  are 
emitted  with  considerable  spatial  variation.    It  can,  therefore,  be 
misleading  to  use  average  emissions  and  concentrations.    It  should  be 
possible  to  develop  "rules",  which  may  be  region-specific,  that  one  would 
expect  to  encounter  local  concentrations  which  are  perhaps  a  factor  of  10 
or  20  higher  than  average. 

These  concentrations  could  be  of  value  in  three  respects.  First,  they 
provide  guidance  about  likely  prevailing  environmental  concentrations  in 
the  various  media  which  result  from  present  or  possible  fctore  emissions. 
This  is  useful  in  designing  analytical  monitoring  and  identificatic:. 
programs . 


Second,  they  provide  guidance  as  to  the  proximity  of  the  concentrations  to 
those  which  are  judged  to  be  of  toxicological  or  aesthetic  significance. 
Generally,  the  estimated  concentrations  will  be  (or  should  be)  well  below 
lethal  concentrations,  but  the  quotient  of  effect-to-estimated 
concentration  provides  a  numerical  expression  of  the  "safety  factor". 
Chemicals  with  low  quotients  are  presumably  of  greatest  concern  and  thus 
require  priority  treatment.  The  combination  of  the  model  with 
toxicological  concentration  information  thus  provides  a  method  of 
establishing  priorities  among  the  large  number  of  chemicals  in  present  use 
and  helps  to  establish  the  relative  priority  of  new  chemicals. 

Third,  the  concentrations  provide  a  starting  point  for  estimation  of  human 
exposure.  Air  inhalation  and  water  ingestion  rates  can  be  combines  with 
the  appropriate  media  concentrations  to  estimate  exposure  by  these  routes. 
The  relative  importance  of  these  routes  then  becomes  obvious. 

The  importance  of  plants  and  vegetation  as  a  source  of  human  and  animal 
intake  of  contaminants  is  now  widely  recognized.  Persistent  organic 
compounds,  such  as  pesticides,  herbicides,  PCEs,  and  dioxins,  have  been 
found  to  partition  into  grains  (31),  fruits  and  vegetables  (26,  32),  pine 
needles  (33),  and  moss  and  lichens  (34).  They  thus  become  directly 
available  to  herbivores  and  humans  through  the  ingestion  of  grains,  fruits 
and  vegetables,  and  indirectly  through  the  consumption  of  domestic  aniT-als 
and  milk. 


Attempts  are  being  made,  notably  by  Travis  and  co-workers  (35),  to  relate 
concentrations  in  plants  and  animals  to  those  of  soil  and  air,  as  estimated 
from  models  such  as  this.  Estimates  can  then  be  made  of  human  intake  by 
ingestion  of  these  foods  in  appropriate  "market-basket"  ratios.  Intake  by 
fish  consumption  can  also  be  included. 

The  exposure  assessment  can  be  extended  further  to  include  a  recently 
developed  steady-state  physiologically  based  pharmacokineetic  model  (36)  to 
predict  partitioning,  accumulation,  and  persistence  in  various  animal  or 
human  tissues  after  prolonged,   continuous  exposure.     Human  tissue 
concentration  data  is  usually  not  available  for  direct  validation,  but 
extrapolation  from  animal  to  human  models  is  possible.  Calabrese  (37)  has 
discussed  the  relevance  of  physiologically-based  pharmacokinetic  animal 
models  for  human  risk  assessment  of  cancer  through  exposure  to  chemical 
carcinogens.     Although   interspecies  and  intraspecies  variation  in 
susceptibility  and  heterogeneity  of  human  response  to  toxic  agents  are 
areas  of  concern,  the  models  can  be  used  to  quantify  target  organ  doses 
between  species. 

Ultimately,  it  is  believed,  a  capability  will  emerge  to  combine 
environmental,  human  exposure,  and  pharmacokinetic  models  in  an  overall 
process  of  tracking  the  pathways  of  a  chemical  including  its  sources,  its 
distribution  in  various  environmental  media,  such  as  air,  water,  sol, 
sediment,  and  food  (as  estimated  here) ,  its  availability  or  exposure  to 
animals  and  humans,  and  finally,  to  physiological  targe*  tissues. 


Finally,  we  believe  that  the  model  can  be,  and  ehould  be,  improved  by 
obtaining  better  estimates  of  certain  transport  and  partltionlno 
parameters.  For  example,  the  mass  transfer  coefficients  and  dlffusivities 
are  average  values  applied  to  all  chemicals.  Oiemlcal  specif  1c  and  region 
specific  values  could  be  used  Instead  to  measure  accuracy.  Host  important, 
the  model  should  be  tested  more  thoroughly  against  a  variety  of  chemicals 
of  widely  differing  properties.  The  principal  difficulties  In 
accomplishing  this  are  the  lack  of  emission  data,  environmental 
concentration  data,  reaction  rate  constants  and  physical  chemical 
propert les . 


CCNCUBICNS 

A  regional  fugacity  model  has  been  developed,  teeted  with  six  chemicals, 
and  found  to  give  a  reasonable  picture  of  environmental  partitioning, 
reaction,  advecticn,  and  transport  characteristics.  Estimated 
concentrations  are  in  order^f-inagnitude  accord  with  reported  values,  but 
there  are  considerable  uncertainties  in  the  emission  rates  on  tfiich  these 
concentrations  depend.  The  estimated  concentrations  maybe  of  value  for 
designing  analysis  programs,  for  setting  priorities,  and  for  asaseseing 
human  exposure  to  present  and  new  chemicals. 
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